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ABSTRACT 

The so-called quasi-Bell entangled coherent states in a thermal environment are studied. In the analysis, we 
assume thermal noise affects only one of the two modes of each state. First the matrix representation of the 
density operators of the quasi-Bell entangled coherent states in a thermal environment is derived. Secondly we 
investigate the entanglement property of one of the quasi-Bell entangled coherent states with thermal noise. 
At that time a lower bound of the entanglement of formation for the state is computed. Thirdly the minimax 
discrimination problem for two cases of the binary set of the quasi-Bell entangled coherent states with thermal 
noise is considered, and the error probabilities of the minimax discrimination for the two cases are computed 
with the help of Helstrom’s algorithm for finding the Bayes optimal error probability of binary states. 
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1. INTRODUCTION 

The choice or design of quantum states for information-processing resource is one of important issues to build 
a useful quantum information system that breaks the performance limit of classical schemes in information- 
communication technologies or that possesses quite new functions. Entanglement may be one of potential 
phenomena and indeed various types of entangled states have been widely discussed from the point of view of 
both theory and experiment and from fundamental properties to its applications. Toward developing macroscopic 
quantum information-communication technologies, we are interested in entangled coherent states.^’^ Above 
all, our interest is focused on the so-called quasi-Bell entangled coherent states in this paper. The quasi-Bell 
entangled coherent states are defined as a set of the following four superpositions of two-mode coherent states: 
|a, ±q;)i 2 -b |-a, TQ:)i 2 and \a, ±a)i 2 - \-a, TQ:)i 2 -^ 

In an ideal situation, it has been clarified that the quasi-Bell entangled coherent states are applicable to several 
applications. Eor example, teleportation and computation by the quasi-Bell entangled coherent states were 
already proposed.The error-free quantum reading scheme was proposed by Hirota,® where phase-shift keying 
(PSK) format was used for its signal modulation. In order to realise such information-communication systems, 
more precise analysis of the states in various realistic situations is needed. In other words, our central concern 
is whether or not the quasi-Bell entangled coherent states are still applicable to information-communication 
technologies even in practical situations. Eor this question, the effect of channel loss to the quasi-Bell entangled 
coherent states was reported by present author in the context of the quantum reading in SPIE Optics-I-Photonics 
2013.® In line with this thought, we aim to clarify the effect of thermal noise for the quasi-Bell entangled coherent 
states. 

This paper is organized as follows. Section 2 reviews the basics of the quasi-Bell entangled coherent states 
when there is no noise. In Section 3, we derive the matrix representation of the density operators of the quasi- 
Bell entangled coherent states in a thermal environment by using the thermalizing operator. The degree of 
entanglement of one of the quasi-Bell entangled coherent states with thermal noise is discussed in Section 4, and 
an optimal quantum discrimination problem for two cases of the binary set of the quasi-Bell entangled coherent 
states with thermal noise is considered in Section 5. In these sections, we execute computer simulations by using 
the matrix representation mentioned in Section 3. Section 6 summarizes our results obtained in this paper. 
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2. QUASI-BELL ENTANGLED COHERENT STATES WITHOUT THERMAL NOISE 

Consider two modes of electromagnetic field with corresponding annihilation operators Oi and 02- Two-mode 
coherent states are defined by \oii,a2)i2 = -Di(ai)-D2(a2)|0,0)i2, where |0, 0)i2 is the two-mode vacuum state 
and Di{ai) = exp[Q;ia| — a*ai] is the displacement operator of the mode i {i = 1 , 2 ). The state |ai, 02)12 can be 
expressed into the form 


101,02)12 = exp[- 


ai| 


|02 


|2 00 00 fi2 


EE 


2 ^ ^ ^ V’^i!«2! 

ni=0n2=0 V -L - 4 , 


|»T-1, '^2)12, 


( 1 ) 


where |ni, 712)12 are two-mode Fock states. The quasi-Bell entangled coherent states are defined by the following 
superpositions of two-mode coherent states.^ 


1^1)12 = Ar+(|/3,/3)i2 + |-/3,-/3)i2), 
1^2)12 = AU(|/3,^)i2-|-/3,-^)i2), 
1^3)12 = A4(|^,-/3)i2 + |-/3,/3)i2), 
1^4)12 = AU(|/ 3 ,-^)i 2 -|-^,^)i 2 ), 


( 2 ) 

( 3 ) 

( 4 ) 

( 5 ) 


where N± = 1 /\/ 2 (l ± exp[— 4 |/ 3 p]) is the normalizing factors. From these expressions of the states, one can 
find the relations 

|4’3)i 2 = (il 0 .R2(7r))|4>l)l2, |^'4)l2 = (il 0 .R2(7r))|$2)l2, (6) 

where R2{tt) = exp[—^0,202] is the 7r-rotation operator. In the context of optical communications, this property 
is understood as the PSK format of laser light. 

The average numbers of photons for the states $1 and ^>3 are given by 


(n(#i)) = ( 7 t.(# 3 )) = 2 |/ 3 ptanh[ 2 |^p]. 


( 7 ) 


It is easy to see that (7i(#i)) = (7 t,(#3)) —>• 0 when /3 — 0 . Similarly, the average numbers of photons for the 
states €>2 and $4 are given by 


(n($2)) = (n($ 4 )) = 2 |/ 3 | 2 coth[ 2 |^| 2 ], (8) 

and (n($2)) = ( 71 (^ 4 )) —t 1 when ,5 —>• 0 . Symbol (n($i)) is abbreviated to ( 71 ) if it is clear from the context. 
The set of the inner products of the quasi-Bell entangled coherent states forms the following Gram matrix. 
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1 
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Ki 3 

0 

I 

0 

12(^4 ^1)12 
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where the entry Ki^ = sech[ 2 |/ 3 p]. From this matrix, it is noticed that the states #1 and $3 are non-orthogonal, 
while any other pair is orthogonal. The orthogonality of the states $2 and $4 is utilized to propose an error-free 
binary PSK quantum reading.® Indeed, the minimum probability of error (in terms of the minimax discrimination 
problem^°“^^) for the states $2 and 4>4 is 


Pg°(4>2,$4) = max 

{P2-I-P4 = 1;P2>0;P4>0} 

while the minimum probability of error for the states and $3 is 

^1 — — 4 pip 3 |i 2 ($l|$ 3 )l 2 


exp[- 4 |/?|^] 

1 -I- exp[- 4 |/ 3 p] ■ 





max 

{pi-l-P3=l;pi>0;p3>0} 


( 11 ) 


















It is well known that the states $2 and #4 are maximally entangled states.^ The degrees of entanglement of 
the quasi-Bell entangled coherent states are given by 


E{\^l)l2) = E{\^3),2) = - 

£'(|$ 2 ) i 2 ) = ^'(1^4)12) = 1 , 


1 + K 


13 


l0g2 
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1 - /^13 
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I-K 13 


( 12 ) 

(13) 


where we have used the entropy of entanglement^^’ defined by 


£^( 1 ^) 12 ) = -Tn {pi log 2 /5i) = -Tr 2 (/32 log 2 P 2 ) (14) 

with the reduced density operators pi = Tr 2 |^*)i 2 (^| and p 2 = Tri|$)i 2 ($|. Eq.(13) shows that ^>2 and ^>4 are 
maximally entangled when there is no noise. 


3. QUASI-BELL ENTANGLED COHERENT STATES IN A THERMAL 

ENVIRONMENT 

In this section, we give the matrix representation of the density operators of the quasi-Bell entangled coherent 
states in a thermal environment. For this aim, we use the thermalizing operator^®’T(0) defined by 

T(0) = exp[—0(aa — a^a^)]. (15) 

In this definition, a is the annihilation operator of the real mode to be primary investigated, and a is the 
annihilation operator of the fictitious mode that satisfy [a, 5,] = 0 = [a, a^] and [a, a^] = 1. Further, cosh0 = 
1/Vl — sinh0 = l/Ve'^A's?’ _ 1 with frequency w, Boltzmann constant fee, and temperature T. By 

using the thermalizing operator, various types of thermal states can be considered: (i) thermal vacuum state, 
|0;O,O) = T(0)|O,O), where |0,0) is the zero-temperature vacuum, (ii) thermal coherent state,|0;ai,Q;2) = 
T( 0 )|ai,a 2 ) = T( 0 )I?(ai)I)(a 2 )|O, 6) (hi) coherent thermal state, \a'i,a'2]9) = .D(Q;j)i)(Qf 2 )T( 0 )|O, 6) (iv) 
thermal squeezed state — thermalizing a squeezed state —^i7-2i squeezed thermal state — squeezing a 

thermal state —gQ ^^ 22,23 thermalizing operator provides a useful tool in the analysis of 

optical quantum states in a thermal environment. 

Here let us describe the model to be investigated in this study. We assume thermal noise affects only mode 
2. Then the quasi-Bell entangled coherent states in a thermal environment is formally written as follows. 


Pl 2 (^l) 

= Tra 

{ (ii 0 t23(0)) |$i)i2|6)3}{i2($i|3(6 

|(il0fj3(0))}' 

(16) 

^ 12 (^ 2 ) 

= Tra 

1 0 T23(^)) ^ 2)12 0)3||i2(^2 3(C 

|(il0tj3(0))}' 

(17) 

^ 12 (^ 3 ) 

= Tra 

1 0 T23(0)^ ^ 3)12 0)3||i2(^3 3(0 

|(il0tj3(0))}' 

(18) 

^ 12 (^ 4 ) 

= Tra 

{ (ii 0 t23(0)) |$4 )i2|6)3}{i2($4|3(6 

|(il0tj3(0))}' 

(19) 


where mode 3 is the fictitious mode for calculating the effect of thermal noise that affects only on mode 2, so 
that we have set T 23 (^) = exp[— 0(0203 — o^Og)]. By using the overcompleteness of two-mode coherent states 
(I/tt^) / d^ 7 id^ 72 | 7 i, 72 )( 7 i, 72 I = I 12 , the four states mentioned above are rewritten as follows. 

= d'7icl^(d272d27^[Ei(/3,/3)F2(/3,/3) +Fi(/3,-/3)F2(/3,-/?) 

+Fi(-/3,/3)F2(-/3,/3) + Fi(-/3,-^)F2(-/3,-/3)]|7i, 72 ) 12 ( 7 ),72l> (20) 


Pl2(^l) 




















Pi2(^2) = ^ 


(/3, P)F2 (/3, /3) - Fi (/3, -/3)i^2 (/?, -/3) 
-i^i(-^,/3)F2(-/3,/3) + i^i(-^,-/3)F2(-^,-/3)ll7i,72)i2(7^,72l, (21) 


Ml 

^ 12 (^ 3 ) = d' 


d27id27;d272dd [Fi(/ 3,/3)F2(-/3, -/3) + Fi(/3, -/3)F2(-/3, /3) 

+Fi(-/3, /3)F2(/3, -/3) + Fi(-/3, -/3)F2(/3, d I 71 , 72 )i 2 ( 7 ( , 72 l> (22) 


^ 12 (^ 4 ) = —X 


d^id^id272d27' [Fi (/3, /3)i^2 (-/3, -/3) - F, (/3, -/?)i^2 (-/3, /3) 

-i^i(-/3,/3)i^2(/?,-/3)+Fi(-/3,-/3)F2(/3,/3)ll7i,72)i2(7(,72l, (23) 


where 


-Fi(/ 3 L,/?ii) = exp 


-|/ 3 p-^l 7 iP+/?L 7 r-^WP + / 3 fl 7 ( 


F2Wl,/3r) = 


-exp 

cosh^ 6 




P*R 


(24) 

72 +72 72 tanh^^' > ( 25 ) 


for Pl^Pr G {—/ 3 ,/ 3 }. From the expressions of Eqs.( 20 )-( 23 ), one can find symmetric relations of the states. 
For Pi2(^i) and ^12(^3), we have ^12(^1) = (li ® .R2(7r))pi2(^3)(ii ® .R2(7r))l. Similarly, ^12(^2) = (li ® 
-R 2 ( 7 r))pi 2 (^ 4 )(ii ( 8 )i? 2 ( 7 r))l for Pi2(^2) and ^12(^4)- As we will see later, this type of symmetry of the states 
reduces an optimization step in the minimax discrimination problem. 

To execute computer simulations, we rewrite the density operators of Eqs.( 20 )-( 23 ) into the corresponding 
matrix representation by Fock states {|n.i, n.2)i2}. The matrix representation of Pi2(^i) ia the basis {\ni, 742)12} 
is given as follows. 

00 00 00 cxD 00 

Pi2(^i) = '^'^[9iW^'rni)9i(.P^ni)923{P,m2,k)9h{^,n2,k) 


mi—0 m 2—0 ni—0 712—0 k—0 


+gi{l3,mi)gl{-P,ni)g23{l3,m2,k)gh{-P:n2,k) 
+gi{-P,mi)gl{/3,ni)g23{-P,m2,k)923iP,n2,k) 
+gi(-/3, mi)gl{-l3, ?4i)g23(-/3, 'rn2,k)9h{-P, ^2, k) 
x|mi, 7772 ) 12 ( 741 , 742 !, 


( 26 ) 


^ 12 (^ 2 ) = 


00 00 CO 00 CXD 

A/)! ^ ^ ^ X! X! [51 (^’"^1)^1 (/?, 7742, %23(^)»T-2,fc) 

mi=0 7712=0 71l=0 712=0 k — 0 

-51 (/3, 774l)g} (-/3, 741)523 (/?, 7742, fc)523(-/3, 742, k ) 
-5 i(-/3 , 7741)5}(/3, 74 i)523 (-/3, 7742, k ) g 23 i ^, 742, k ) 
+5i(-/3, 744i)5}(-/3, 74i)523(-/3, 7742, k ) g 23 {- P , 742, k ) 
XjTTll, 7742)12(741,742!, 


( 27 ) 

















Pi 2(^3) = ^ ^ ^ ^ ^[5l(/3,"il)ffr(^,?T-l)523(-/3,m2,%23(-/3:'«2,fc) 

mi —0 7712=0 ni =0 712=0 fc =0 

+ffi(/3, rni)gl{-P, ?T-i)ff23(-/3, ^2, k)g23{P, n2,k) 

+9i{-P, 'mi)gl{P, rii)g23iP, ^2, fe)ff23(-/3) ^2, fc) 
mi)5i (-^, ?T-i)ff23(/3, TO2, k)g23{l3, n2,k) 
x|mi,m2)i2(ni,n2|, (28) 


Pi2(^4) = X! X! X! X! X![5l(/5’'^l)5t(/3,'rai)523(-/?,?T42,fc)g23(-/?)»T-2,fc) 

mi=0 7712=0 711=0 712=0 fe = 0 

-ffi(/3, mi)5i (-/?, ni)g23{-P, m 2 , fc)ff23(/3> n2,k) 

-gi{-P, mi)gl{P, 711 ) 323 (/ 3 , m 2 ,k)g 23 {-P, n2,k) 

+3i(-^, mi)gl{-(3, 711 ) 323 (/?, ^2, k)g23{l3, n2,k) 

x|777i, 7712 ) 12 ( 771 , 712 !, 


where 


5i(/3',fei) 

923{P', ^ 2 ,^ 3 ) 


10/12 o/fci 

i(fci|/ 3 ')i=exp[-^]-^ 


2 

23(^2, fc3|T23(0)|/3^ 6)23 
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cosh0 




0 , if k2 < ks 


(29) 


(30) 


(31) 


and we have applied the calculation result of Ref. 24 to obtain Eq.(31). From Eqs.(26)-(29), the photon-number 
distributions P{ni,n 2 ) = (tii, 772|Pi2(^i)|77i, 712 ) can be derived. Figure 1 shows examples of photon-number 
distribution P{ni,n 2 ) of the state $2 for /? = 2.0. The graphs (a-1) and (a-2) show the case of 0 = 0, that is, 
noiseless case, and the graphs (b-1) and (b-2) show the case of 0 = 0.5. When there is no thermal noise (0 = 0), 
the distribution P{ni, 772 ) is not equal to zero only if 771 -I -772 is odd. From this fact, |$ 2 )i 2 is understood as the 
odd two-mode coherent state.On the other hand, we see that this property is destroyed due to the effect of 
thermal noise from the case of 0 = 0.5, in which the dispersion of the photon-number distribution is appeared in 
the direction of the 772 -axis. 


4. EVALUATION OF THE DEGREE OF ENTANGLEMENT OF $2 

As mentioned in Section 2, the states $2 and $4 are maximally entangled when there is no noise (Eq.(13)). 
Further, these two states are associated with each other by a unitary operator (Eq.( 6 )). In this section, we 
consider the entanglement property of the state $2 with thermal noise. The analysis method of the entanglement 
property of ^ 12 (^ 2 ) is based on Ref. 26. 

The fully entangled fraction for the state d‘i 2 is defined as 


/(d-12) = maxi2(e|d-i2|e)i2, 

where the maximum is taken over all maximally entangled states |e)i2. Here we let 


/(Pi 2(^2)) =max i2(^2(a)|Pi2(^2)|^2(a))i2 


(32) 


(33) 


in order to analyse the state p;^ 2 (^ 2 )- Note that we chose the state |$ 2 (q!))i 2 = M-{a) [|a,a)i 2 — |—a, — 0 ) 12 ] 
as a maximally entangled state, where Af-{a) = 1/ 1 / 2(1 — exp[—4|ap]). The complete form of Eq.(33) is shown 












in Appendix A. By using the resulting value of f {^ 12 (^ 2 )), a lower bound of the entanglement of formation is 
given by 

E{p,2i^2)) > H(/(Pi2 (^2))), (34) 

where 

with the binary entropy function h2{x) = —xlog 2 x — {1 — x) log 2 (l — x). 

Figure 2 shows the numerical behavior of f {^ 12 (^ 2 )) and the lower bound of ^(^,^ 2 (^ 2 )) obtained from 
the value of f{pi 2 {^ 2 ))- Although it might not be tight, the degradation of the degree of entanglement due 
to thermal noise is obviously observed through the graph (2) of Figure 2. Thus the degree of entanglement 
is degraded by increasing the amplitude parameter /3 and by increasing the thermal noise parameter 6. The 
threshold level of / in the function H(/) is 0.5. Hence, from the graph (1) of Figure 2, we roughly estimate that 
the limit of the thermal noise parameter 0 to keep entanglement is ^ 0.7. 
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Figure 1. Photon-number distribution P{ni,n 2 ) of the state €> 2 . (a-1) P(ni,n 2 ) for 6 = 0 and /? = 2.0. (a-2) Contour 
plot of (a-1). (b-1) P{ni,n 2 ) for 9 = 0.5 and /3 = 2.0. (b-2) Contour plot of (b-1). 
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Figure 2. Entanglement property. (1) behavior of /(Pi 2 (^ 2 )). (2) lower bound of ^(^^ 2 (^ 2 )) obtained from f{pi 2 {^ 2 ))- 


5. ERROR PERFORMANCE IN THE MINIMAX DISCRIMINATION PROBLEM 

Suppose there are two quantum states ctq and di, where di > 0 and Trd^ = 1. If the probability distribution 
P = {Po,Pi) of the states is given in advance, one can employ the quantum Bayes strategy.At that time, 
the average probability of error of the Bayes optimal discrimination for given p is formally written as 


pmin/ ) ^ 

n 


PoTrlTido+PiTrI7o^i , 


(36) 


where U = {7Jo,/7i : IJq + 77i = 1 and Ili > 0} is a positive operator-valued measure (POVM) which describes 
a detection profile for the states. 

Clearly .Pe™™(p) depends on the choice of p. To remove the dependency of choice of the probability distribution 
p, we employ the quantum minimax strategy.When the quantum minimax strategy is employed, the average 
probability of error in the minimax optimal discrimination is given by 


P° = min max 
n p 


poTrlZido -|-piTr77odi 


= max mm 
P 77 t 


poTr/Zido + piTrlloai 


= max [^ 3 “” 


(P)] ■ (37) 


It is not easy to find the optimal probability distribution p° in general. However, if the states do and di are 
symmetric, the optimal distribution becomes the uniform distribution u = ( 1 / 2 , 1 / 2 ). At that time, we have 
P° = i/f""(l/2,1/2). Recalling the Helstrom’s algorithm for calculating the minimum error probability for 
binary states,it becomes 

^ ^ E (38) 

Ai>0 

where Xi are the eigenvalues of do — di. 

Applying the result of Eq.(38), we can compute the error probabilities for the set of the states $2 and €> 4 , and 
for the set of and $3 in terms of the minimax discrimination problem. Results of the numerical calculation 


















are shown in Figure 3. The graphs (a-1) and (a-2) illustrate and Pg°(0; # 3 ), respectively, 

for 9 = 0.01,0.05,0.1,0.3,0.5,0.7, and 0.9. Comparing these two graphs, it is noticed that the behavior of 
.Pg°(d; $ 2 , ^ 4 ) is different from that of .Pg°(0; ^> 1 , $ 3 ). Recall that .Pg°(0; $ 2 , ^ 4 ) = 0 if d = 0. As a value of 
the thermal noise parameter 9 increases, the error probability Pg°( 0 ;$ 2 ,^ 4 ) rapidly moves upward, while the 
error probability .Pg°( 0 ; $ 1 , $ 3 ) changes moderately as if to be the case of binary single-mode coherent states 
{|/3)i, |—/3)i}. The graphs (b-1) to (b-6) involve numerical comparison between P°{9-, $ 2 , ^ 4 ) and Pg{9- ^> 1 , ^> 3 ) 
for each 9. In all cases from (b-1) to (b- 6 ), Pg°( 0 ;# 2 ,^ 4 ) is always smaller than Pg°(0; ^> 1 , $ 3 ). But, the gap 
between P°{9-, $ 2 , ^ 4 ) and P°{9; $ 1 , # 3 ) is getting smaller when the thermal noise parameter 9 increases. When 
9 is small enough {9 ~ 0.1 or less), one can easily find a clear gap between P°{9-, # 2 , ^ 4 ) and P°{9-, $ 1 , $ 3 ) in 
the region 0 ~ (n) ~ 8. When 9 increases, Pe{9-, # 2 , ^ 4 ) approaches to Pe{9; $ 1 , $ 3 ), and finally Pg{9] ^> 2 , ^> 4 ) 
almost touches to Pg{9] # 1 , # 3 ). 

5.1 Conclusions 

The so-called quasi-Bell entangled coherent states, # 1 , $ 2 , ^3 and $ 4 , were studied in the presence of thermal 
noise. In the analysis, we assumed thermal noise affects only one of the two modes of each state. By using 
the thermalizing operator with thermal noise parameter 0 , the matrix representation of the density operators 
of the quasi-Bell entangled coherent states with thermal noise was derived. To seek the entanglement property 
of the states in the presence of thermal noise, a lower bound of the entanglement of formation for the state 
$2 was numerically computed. From this, we observed that the range of 9 to keep entanglement of $2 is 
limited to a small region. We also considered the optimal quantum discrimination problem for the states ^2 
and $4 and for and $ 3 . From numerical comparison, it is noticed that the behavior of the error probability 
Pg°(0;$2,^4) is different from that of .Pg°(0; $ 1 , # 3 ), and the gap between P°{9;^2,^4) and .Pg°(0; $ 1 , ^> 3 ) is 
getting smaller when 9 increases. The advantage of the use of the states of $2 and 4>4 against the states of 
and ^>3 vanishes when 9 is large. In summary, the effective range for the use of the states of $2 and ^>4 in 
information-communication applications would be limited rather than other two, so that careful design for such 
applications is needed. 
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Eq.(33) is given by 
/(Pi 2 (^ 2 )) = max 
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This maximization problem with respect to a was numerically solved to draw the graph (1) of Figure 2 under 
the restriction that all parameter is real. 



























